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Abstract We present an initial look at the far infrared-radio correlation within the star- 
forming disks of four nearby, nearly face-on galaxies (NGC 2403, NGC 3031, 
NGC 5194, and NGC 6946). Using Spitzer MIPS imaging and WSRT radio 
continuum data, we are able to probe variations in the logarithmic 70 /im/22 cm 
(g7o) flux density ratios across each disk at sub-kpc scales. We find general 
trends of decreasing q7o with declining surface brightness and with increas- 
ing radius. We also find that the dispersion in q-jo within galaxies is compa- 
rable to what is measured globally among galaxies at around 0.2 dex. We have 
also performed preliminary phenomenological modeling of cosmic ray electron 
(CRe~) diffusion using an image-smearing technique, and find that smoothing 
the infrared maps improves their correlation with the radio maps. The best fit 
smoothing kernels for the two less active star-forming galaxies (NGC 2403 and 
NGC 3031) have much larger scale-lengths than that of the more active star- 
forming galaxies (NGC 5194 and NGC 6946). This difference may be due to 
the relative deficit of recent CRe~ injection into the interstellar medium (ISM) 
for the galaxies having largely quiescent disks. A more complete discussion of 
this proceedings article can be found in Murphy et al. (2005). 
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1. Introduction 

Soon after the Westerbork Radio Synthesis Telescope (WSRT) went online, 
van der Kruit (1971) made the first report of a correlation between 1415 MHz 
radio and 10 /im mid-infrared luminosities for a sample of Seyfert galaxy nu- 
clei. This discovery was soon extended to the nuclei of normal spirals, though 
the dispersion was found to be smaller in the correlation for the original sam- 
ple of Seyferts (van der Kruit (1973)). It was not until the coming of IRAS 
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that the optically thin radio continuum emission from galaxies was found to 
be better correlated with the far infrared (FIR) dust emission of galaxies with- 
out an active galactic nucleus (AGN) (de Jong et al. (1985); Helou, Soifer, & 
Rowan-Robinson (1985)). 

The connection between radio and infrared emission from galaxies is that 
they are both powered by massive stars, as pointed out originally for starbursts 
by Harwit & Pacini (1975). Young massive stars, which heat up dust to provide 
the FIR emission, are the same stars which end as supernovae (SNe) and bring 
about the synchrotron emission observed at radio wavelengths. If this general 
picture is correct, the fact that the mean free path of UV photons (~100 pc) 
which heat the dust is much less than the diffusion length for a CRe~ (~1- 
2 kpc) suggests that the radio image should resemble a smeared version of 
the infrared image. This idea was first introduced by Bicay & Helou (1990), 
who attempted to model the propagation of CRe~s by smearing IRAS data 
of galaxies using parameterized kernels containing the physics of the CRe 
propagation and diffusion, to better match the morphology of the correspond- 
ing radio data. 

In an attempt to better understand the FIR-radio correlation, we are using 
MIPS infrared data from the Spitzer Infrared Nearby Galaxies Survey (SINGS) 
legacy science project (Kennicutt et al. (2003)). In this initial study, we exam- 
ine the variations of the FIR-radio correlation on sub-kpc scales within four 
of the nearest face-on galaxies in the SINGS sample: NGC 2403, NGC 3031 
(M81), NGC 5194 (M51a), and NGC 6946. Using high resolution Spitzer 
imaging, we are also able to test the smearing model of Bicay & Helou (1990) 
with greater accuracy, at higher spatial resolution, and in more galaxies, with 
the aim of gaining better insight into CRe diffusion and confinement within 
galaxy disks. For a complete discussion, we direct the reader to Murphy et al. 
(2005). 

2. Results 
q 70 Maps 

In Figure 1 we plot q 70 = log ( (22 cm) [jy] ) ^ or eacn § a l axv - We find 
that elevated q-jo ratios are associated with bright structures in the infrared and 
radio images of each galaxy. The most obvious cases are seen for the bright 
spiral arms of NGC 3031, NGC 5194 and NGC 6946. All three galaxies show 
enhanced qjo ratios along their arms, with local peaks centered on H II regions, 
and depressed ratios located in the quiescent inter-arm and outer-disk regions 
of each galaxy. For NGC 2403, which does not have a grand-design spiral 
morphology, we still see q-j® peaks associated with H II regions. 

To quantify radial variations in q-jQ, along with any dependencies of qjo on 
the 70 ^m surface brightness within each galaxy, we perform aperture pho- 
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Figure 1. 570 maps for pixels having a 3-er detections in both the input radio and 70 
maps. The q^o ratios are enhanced along spiral arms and on H II regions due to 70/im emission 
being more strongly peaked than the 22 cm emission. 



tometry using apertures having diameters equal to the FWHM of the MIPS 
70 /im PSF (i.e. 17"). The measured dispersion of qjQ within each galaxy 
is less than ~0.25 dex. We find a general trend of increasing qjo ratios with 
increasing 70 /*m surface brightness and decreasing radius. The residual dis- 
persion around the trend of (770 with increasing radius is found to be larger than 
the residual dispersion around the trend of qjo with increasing 70 ^m surface 
brightness by ~0. 1 dex. This suggests that the star formation sites within the 
disk are more important in determining the overall appearance of the qjo maps 
compared to the exponential profiles of the disks themselves. 

Infrared/Radio Relations Inside and Among Galaxies 

In Figure 2 we show our local sub-kpc qjo ratios together with global qjo 
values for 1752 galaxies from the study of Yun, Reddy, & Condon (2001). 
Using their cataloged 1.4 GHz NRAO VLA Sky Survey (NVSS) and 60 and 
100 fim IRAS data, we converted IRAS based q^o values to Spitzer 570 values 
using the SED models of Dale & Helou (2002). We find that the dispersion in 
the local sub-kpc qjo ratios within our sample galaxies is comparable to that 
of global ratios. In the Yun, Reddy, & Condon (2001) sample, qjo appears to 
be roughly constant with increasing galaxy luminosity while, within each disk, 
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qjo ratios clearly increase with the 70 fim luminosity. This difference in the 
observed trend between q^o and luminosity within and among galaxies can be 
explained by the diffusion of CRe~ within the galaxy disks. 
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Figure 2. Matched 1.5 kpc aperture measured 570 ratios for each sample galaxy plotted with 
global Q70 ratios estimated from the data presented by Yun, Reddy, & Condon (2001). While 
the global qro ratios are constant with 70 ^im luminosity, the local sub-kpc 570 ratios are not. 



3. Image-Smearing and Cosmic Ray Diffusion 

To help understand the underlying physics of the correlation, we look to 
see whether the image-smearing model works to improve the correlation be- 
tween the infrared and radio morphologies of galaxies. Determining the func- 
tional form of the best fit smearing kernel provides insight into the propagation 
and diffusion characteristics of cosmic rays within galaxy disks. In this study 
we examine the behavior for Gaussian and exponential kernels oriented in the 
planes of the sky and each galaxy disk. 

We define the residuals between the smeared 70 /xm map ( I(t, p, l):t =shape, 
p =orientation, and Z=scale-length) and the observed 22 cm map (R) by, 

<f>(Q,t,p,l) = , (Ij 

where Q = and j indexes each pixel. In Figure 3 we plot log ((/>), 

as a function of the scale-length (/), for each galaxy and find that the image- 
smearing technique improves the overall correlation between the radio map 
and the 70 /im image by an average of ~0.2 dex. Unlike Marsh & Helou 
(1998), we find exponential kernels are preferred independent of projection 
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Figure 3. Residuals of observed radio maps with smeared 70 /im images, denned by 
4>(Q, t, p, I), as a function of smearing scale-length shown for each kernel per galaxy. Residuals 
for the active star-forming galaxies are minimized by kernels having < 1 kpc scale-lengths while 
the less active star-forming galaxies are better fit by those with scale-lengths > 1 kpc. 



over Gaussian kernels, suggesting that additional processes such as escape and 
decay appear necessary to describe their evolution through the galaxy disks. 

We find distinct behavior between our active and less active star-forming 
galaxies. Figure 3 shows that scale-lengths <1 kpc work best to improve the 
correlation for the more active star-forming galaxies (NGC 5194 and NGC 6946: 
SFR>4 M Q /yr), while scale-lengths >1 kpc work best for the less active star- 
forming galaxies (NGC 2403 and NGC 3031: SFR <lM©/yr). We also find 
that by smearing the 70 /um images, we are able to reduce the slope of the ob- 
served non-linearity in qjo with 70 /im surface brightness by >25% for each 
galaxy. This suggests that the non-linearity may be due to the diffusion of 
CRe~s and is not expected to be found in the global correlation when integrat- 
ing the flux over entire galaxies. 



4. 



Conclusion 



In examining the behavior of q-jo within the disks of our sample galaxies, 
along with our image-smearing analysis, we find: 
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1 An empirical trend of q-jQ decreasing with declining 70 fim surface bright- 
ness and increasing radius to be a general property within the galax- 
ies. However, the dispersion measured in q-jQ at constant 70 fim surface 
brightness is found to be smaller than at constant radius by ~0.1 dex, 
suggesting that star formation sites are more important in determining 
the qjo disk appearance than the underlying exponential disks. 

2 The dispersion in the global FIR-radio correlation is comparable to the 
dispersion for q-jQ on sub-kpc scales within the galaxy disks. Also, the 
trend of increasing q-?o ratio with increasing 70 fxm luminosity within 
each galaxy is not observed in the global correlation. 

3 The phenomenological modeling of cosmic ray electron (CRe~) dif- 
fusion using an image-smearing technique is successful as it both de- 
creases the measured dispersion in and reduces the slope of the ob- 
served non-linearity in qjo with 70 surface brightness. The latter 
suggests that the non-linearity may be due to the diffusion of CRe~s 
from star-forming regions. 

4 Exponential kernels work marginally better to tighten the correlation 
than Gaussian kernels, independent of projection. This result suggests 
that CRe evolution is not well described by a pure random-walk diffu- 
sion and requires additional processes such as escape and decay. 

5 Our two less active star-forming galaxies require kernels having larger 
scale-lengths to improve the correlation while our two more active galax- 
ies require smaller scale-lengths. This difference may be due to time- 
scale effects in which there has been a deficit of recent CRe injection 
into the ISM of the two less active star-forming galaxies, thus leaving 
the underlying diffuse disk as the dominant structure in the morphology. 
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